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SECTION  I 


INTRODUCTION 


1.  Increased  AF  ur;a  of  hydraxina  propellants  and  more  stringent 
Occupational  Safety  and  Health  Administration  (OBHA)  and  Environmental 
I’roteotion  Agency  (ERA)  regulations  have  greatly  increased  the  need 

to  pradiat  concentrations  of  polluting  chemicals  discharged  into  the 
environment.  In  January  1976  the  AF  Space  and  Missile  Systems  Organisa- 
tion (SAMSO)  identified  Technology  Need  (TN-SAMSO-CEEDO-2109-76-45) 
that  addressed  this  issue.  The  TN  opeaifioally  stated  that  "The  fata 
of  the  hydraxino  fuels  entering  the  environment  is  insufficiently 
characterised"  (Reference  1) . Prediction  techniques  available  at  missile 
nitns  prior  to  this  study  were  not  easily  adaptable  because  they  re- 
quired temperature  difference  input  data  not  commonly  available  and 
beaause  photochemical  reactions  and  atmospheric  decompoeition  were 
not  considered. 

2.  Hydrasme  propellants  are  ueed  in  a variety  of  systems , such  as 
■pact  launch  vehicles,  strategic  missiles  and  astellitss  and  the  F-16 
Emergency  power  Unit.  The  environmental  impact  of  hassrdous  chemicals 
vented  into  the  atmosphere  or  accidentally  spilled  during  loading  and 
off-loading  operations  oan  be  more  quickly  elucidated  through  use  of 
computerised  evaporation  end  dispersion  analysis  tsohniquss.  Thsss 
predictive  techniques  ere  crucial  to  asseasing  the  hasard  hydraainas  may 
posn  to  man. 


3.  The  objective  of  this  technical  report  is  to  present  a taohniqus 
for  computing  evaporation  rates  from  accidental  spills  and  to  illuatrata 
how  « simple  physical  dispersion  modal  can  ba  used  to  maka  air  quality 
predictions  downwind  from  the  spill.  This  report  assumes  the  evaporation 
from  a non-porous,  non-absorbing,  flat  surface  and  does  not  address  atmos- 
pheric ahemical  reactions  of  the  evaporating  propellant.  In  addition,  the 
evaporation  model  pertains  strictly  to  purs  missile  propellants  end  cannot 
be  applied  straight  forwardly  to  mixtures,  such  as  used  in  TITAN  missiles 
end  the  F-16  aircraft.  The  model  is  currently  being  modified  to  handle 
propellant  mixtures. 

The  computation  of  evaporation  rates  and  downwind  concentrations 
are  important  sinas  they  have  a major  bearing  on  the  else  of  the  hesard 
corridor  that  must  be  evacuated  to  protect  the  populace  from  the  epill. 
Example  hasard  corridor  computations  for  drum,  trailer  and  railcar  pro- 
pellant spills  are  illustrated  in  Section  XV. 

4.  A number  of  dispersion  models  were  reviewed  to  determine  which  wee 
most  applicable  to  the  propellant  problem.  They  include  the  box  model 
commonly  used  for  urban  area  aourcas,  the  Gaussian  plume  model  widely 
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used  for  continuous  elevated  and  ground  sources , the  Gaussian  puff 
model  used  for  instantaneous  explosive  sources  and  the  empirical  tempera- 
ture difference  model  used  at  the  Kennedy  space  Center  and  Vandanbarg 
APB,  California  for  safety  prediction  at  missile-test  ranges.  The  con- 
ventional Gausaian  model  was  selected  for  this  report.  The  rationale 
for  this  selection  is  presented  in  Section  IV,  which  also  inoludes  dis- 
cussion of  the  calculation  technique  for  determining  downwind,  ground- 
level  concentrations  and  the  hasard  corridor  to  be  evaouated. 
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SECTION  II 


DISCUSSION  OF  SOURCES 


Discharge  of  the  three  hydrasine  missile  fuel*,  anhydrous  hydraaina 
(n2H4> , uneymmetriael  dimethylhydraaine  (UDMH)  and  monomethylhydraeina 
(MMH) , Into  the  environment  can  occur  at  Air  Foroa  fuel  storage  and 
transfer  facilities,  TITAN  XI  strategic  missile  sites  or  the  TITAN  III 
spaae  launch  vehicle  operational  areas  at  the  Eastern  and  Western  Test 
Ranges.  In  addition,  hydraaina  and  monomethylhydrasine  spills  are  a 
potential  problem  with  future  apace  transportation  systems,  aircraft 
auxiliary  power  systems  and  starter  cartridges  that  utilise  hydraaina 
propellants.  About  2.4  million  kilograms  (5.2  million  pounds)  of 
propellants  are  moved  over  320,000  km  (200,000  miles)  annually 
(Reference  2) . In  addition,  large  quantities  are  Btored  at  facilities 
such  as  Rocky  Mountain  Arsenal,  Edwards  AFB  and  Aerojet  General  Corp, 

CA.  Small  quantities  are  also  stored  at  user  bases  all  over  the  country. 
Shipments  are  made  in  trailers,  railcars  and  drums.  The  quantities 
of  chemicals  and  the  various  shipping  containers  are  listed  in  Table  1 
(Reference  2) . 

It  is  estimated  that  0.6  million  kg  of  N2H4,  0.23  million  kg  of  MMH 
and  0.64  million  kg  of  UDMH  are  moved  annually  within  the  United  States 
(Reference  2) , The  potential  for  a spill  during  movement  of  propellant 
from  the  manufacturer  to  the  storage  facility  and  finally  to  the  user 
therefore  exists. 


TABLE  1.  PROPELLANT  SHIPPING  CONTAINERS  AND  QUANTITIES* 

Liters 


Propellant 

Trailer 

Rail  Car 

Drum 

n2«4 

18000 

22,000 

200 

MMH 

21000 

36,000 

190 

UDMH 

21000 

36,000 

180 

‘Reference  2 
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SECTION  til 


PROPELLANT  EVAPORATION 
1.  Computational  Procedure 

To  utilize  atmospheric  dispersion  and  chemical  rate  equations 
for  the  calculation  of  pollutant  concentrations  in  the  air,  the  amount 
of  pollutant  released  into  the  atmosphere  per  unit  time  must  first  be 
quantified.  The  evaporation  rate  of  a propellant  accidentally  spilled 
is  a function  of  ambient  air  temperature,  wind  speed,  solar  radiation, 
dimensions  of  the  propellant  spill  and  its  volatility  and  diffusion 
characteristics  (Reference  3) . In  this  analysis,  three  assumptions 
are  made  concerning  the  evaporation  ratei 

a.  The  spill  process  is  adiabatic  and  chemically  stable. 

b.  The  propellant  is  spilled  as  a liquid  without  atomisa- 
tion. 

a.  Evaporation  occurs  at  a steady-state  pool  temperature. 


The  evaporation  rate  was  oomputad  by  the  method  of  Mackay  (Refer- 
ence 3)  and  is  a function  of  concentration  driving  force,  as  determined 
from  vapor  pressure)  mass  transfer  rate,  as  determined  by  wind  generated 
turbulence)  and  diffusive  properties  of  the  air-liquid  intarfaoa.  The 
methodology  assumes  evaporation  of  a pure  liquid  and  ideal  gaa  behavior 
of  the  film. 

The  rate  of  evaporation  of  a liquid  pool  can  be  described  by  the 
following  equation) 


qm  " VV(mw)Vrtp  (1) 

where 

Qw  - mass  transfer  rate  of  propellant  into  the  air,  kg/hour. 

km  ■ mass  transfer  coefficient,  m/hr. 

Pv  - vapor  pressure  of  propellant,  kpa. 

Tp  x equilibrium  pool  temperature,  °K. 

R x universal  gas  constant,  8.314  kPa^m^/kg’moie'OK. 

MW  x molecular  weight  of  propellantB 

A ■ Area  of  spill,  m2. 

P 
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The  mass  transfer  coefficient  (k  ) is  a function  of  the  mean  wind 
speed  (0) , a constant  (n)  related  to  ground  roughness  end  the  tempera- 
ture profile  of  the  atmosphere,  and  finally  the  Schmidt  number  defined 
in  Appendix  C (Reference  3) . Although  wind  speed  can  be  uasily  measured 
and  the  roughness  factor  assumed  for  average  atmospheric  condictions,  the 
Schmidt  number  must  be  specifically  computed  for  the  propellant-air  ayatem 
at  the  equilibrium  pool  temperature. 

The  pool  temperature  (T  ) is  determined  from  the  steady  state  energy 
balance  equations 

Qq  + + Ss  * Qa  - Qe  + he  <2) 

whoru 

Qg  » the  convective  heat  transfer  from  the  ground  and  pool,  J/hr 

Qh  ■ the  convective  heat  transfer  from  the  atmosphere,  J/hr 

Qs  - solar  inaolation,  J/hr 

Qa  - radiative  heat  transfer  from  the  atmosphere,  J/hr 

Qe  - radiative  heat  amission  from  fchs  pool  of  liquid,  J/hr 

H£  - heat  loss  from  the  pool  due  to  evaporative  cooling,  J/hr 

The  basic  assumption  regarding  evaporation  is  that  the  liquid  pool  tempera- 
ture ia  in  a quasi  steady-state.  The  calculated  pool  temperature  is 
defined  as  that  equilibrium  temperature  which  results  in  a zero  energy 
exchange  between  the  pool  and  the  environment.  The  pool  temperature 
that  satiafiee  the  energy  balance  must  be  computed  by  trial  and  error. 

Aa  a atarting  point  for  the  computations,  the  liquid  pool  tempera- 
ture ia  assumed  equal  to  the  air  temperature.  Baaed  upon  this  assumption, 
the  sir-propellant  vapor  film,  the  liquid  propellant  properties  and  the 
mass  transfer  rate  coefficient  are  calculated.  The  Newton-Raphson  itera- 
tive method  is  used  to  calculate  the  specific  pool  temperature  which 
satisfies  the  energy  balance,  iterative  computations  continue  until  the 
change  in  successive  pool  temperatures  is  less  than  0.001°c.  The  pool 
temperature  will  rarely  be  equal  to  the  air  temperature  due  to  radiative 
heat  transfer  and  evaporative  cooling.  Mackey  states  that  "Clearly,  it 
is  inacourata  to  assume  that  the  pool  temperature  will  equal  the  air 
temperature"  (Reference  3) . 

The  pool  temperature  computed  by  the  Newton-Raphson  msthod  is  thsn 
compared  to  the  originally  assumed  pool  temperature  upon  whioh  the  vapor 
and  liquid  pool  propartiea  and  mass  transfer  rats  coefficient  were  cal- 
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culated.  If  the  difference  is  greater  than  0.1°C  a fractionally  modified 
temperature  is  calculated  in  a second  iteration  set  and  new  vapor  and 
liquid  propellant  properties  and  a new  new  mass  transfer  rate  coefficient 
are  calculated.  These  are  then  used  for  another  solution  to  the  equili- 
brium pool  temperature.  The  procedure  is  continued  until  successive 
values  differ  by  less  than  0,1°C  (See  the  aomputer  flow  chart  in 
Appendix  B) . A detailed  description  of  the  pool  temperature  and  evapora- 
tion rate  equations  is  presented  in  Appendix  A. 

2 . Evaporation  Bates 

Figures  1 through  3 illustrate  NjH,;,  MMH  and  UDMH  evaporation  rates 
for  typical  trailer,  railcar  or  drum  transportation  conditions.  Dashed 
lines  specify  rates  used  to  illustrate  sample  dispersion  problems  in 
Section  TV.  Evaporation  rates  were  calculated  for  three  ambient  air  tem- 
peratures using  the  computer  program  described  in  Appendices  A and  B with 
the  following  acsuiqptionsi 

s.  The  ground  temperature  ie  equal  to  the  ambient  air  temperature. 

b.  wind  speed  is  assumed  constant. 

c.  The  daily  maximum  solar  insolation  (R  ) rates  used  are  2.5 (IQ6) 
J/m2-hr  for  a 0°C  clear  winter  day,  3.1(10®)  J/m2*hr  for  a 15°C  clear 
spring  day,  3.8(106)  j/m2.hr  for  a 30°C  claar  summer  day. 

d.  The  depth  of  spill  is  assumed  constant  at  2.5  om. 

3.  sanaitivity  Analysis  of  Evaporation  Input  Data 

Nina  input  parameters  wars  selectively  varied  to  assess  their  effects 
on  ths  propellent  evaporation  rate.  Table  2 summarises  the  parameters 
and  presents  ■ qualitative  evaluation  o£  their  effect  on  the  propellant 
evaporation  rate.  The  model  is  most  sensitive  to  spill  area,  ground 
temperature,  ground  roughness  and  propellant  type.  Solar  insolation, 
wind  speed  and  air  temperature  represent  medium  sensitive  parameters  as 
does  ground  temperature  when  it  is  below  the  pool  temperature  and  heat 
transfer  by  conduction  occurs.  When  solar  insolation  is  below 
0.8  MJ/m2*hr,  typical  of  times  one  to  two  hours  after  sunrise  end  one  to 
two  hours  before  sunset,  evaporation  is  insensitive  to  insolation.  Vari- 
ations in  pool  depth  and  atmospheric  emissivity  cause  no  significant 
change  in  the  evaporation  rate. 

Senaitivy  analyses  indicate  that  an  error  in  spill  area  results  in  a 
proportional  error  in  the  evaporation  rate.  Since  irregular  spill  areas 
may  be  difficult  to  measure,  assessment  of  this  parameter  may  result  in 
ths  largest  error  to  ths  evaporation  rate.  A 10°C  error  in  the  ground 
tenperaturc  (T0)  results  in  a 24  percent  error  in  the  evaporation  rate 
when  T0  is  below  the  pool  temperature.  When  Tq  exceeds  ths  pool 
temperature,  10°c  error  in  T^results  in  a 55  percent  error  in  the  evapor- 
ation rate.  A variation  of  the  roughness  factor  (n  is  assumed  0.25  for 
turbulent  atmospheric  conditions,  AT  « -4 . 18°c/1000  feet)  between  611  end 
0.5  results  in  evaporation  ratee  between  1220  and  140  kg/hour, 
respectively.  Observed  values  of  n have  ranged  from  0,04  <n<0.93. 


6 


Figure  2.  MMH  Propellant  Evaporation  Ratee  ae  a Function  of 
Ambient  Air  Temperature  and  Spill  Volume. 


SMIL  VOLUMI-UrfftS 

Figura  3.  UDMH  Propallant  Evaporation  Rataa  as  a Function  of 
Ambiant  Air  Tamparatura  and  Spill  Voluma. 
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Analyses  indicate  that  evaporation  rates  can  be  an  order  of 
magnitude  different  from  those  obtained  using  the  same  evaporation 
model  without  heat  gain  and  losses  to  the  evaporating  pool.  For 
example,  a hydrazine  spill  covering  an  area  of  74S  m«  evaporates  at 
615  kg/hr,  rather  than  170  kg/hr,  when  heat  load  input  data  are  con- 
sidered. Thia  example  was  calculated  for  a windspead  of  2.6  m/s, 

Tg  ■ T»  ■ 17.8°C  and  Rg  - 3.8  MJ/m2-hr.  In  general,  one  may  state 
that  the  evaporation  rates  of  UDMH  and  MMH,  for  similar  meteorologi- 
cal and  spill  conditions,  are  eleven  and  three  times  that  of 
hydrazine,  respectively.  The  detailed  sensitivity  analysis  is  pre- 
sented in  Appendix  C. 

Comparison  of  predicted  and  observed  evaporation  rates  show 
good  agreement  for  initial  (first  hour  average)  evaporation  rates. 
Experiments  were  conducted  in  fume  exhaust  hoods  under  controlled 
air  flows  and  room  temperatures  and  outdoors,  during  both  cloudy 
and  sunny  days.  Results  indicate  that  the  model  fails  to  predict  an 
observed  decrease  in  evaporation  rata  with  time.  The  discrepancy 
appears  to  be  due  to  the  absorption  of  and reaction  with  atmospheric 
carbon  dioxide  and  water  vapor  (Reference  4) 


TABLE  2.  QUALITATIVE  SENSITIVITY  EVALUATION  OP 
EVAPORATION  RATE  TO  INPUT  VARIABLES 


Input  Variables  Sensitivity 


Symbol 

Namo 

High 

Medium 

Low 

KP* 

Propellant  Type 

+ 

tg 

Ground  Temperature 

+ 

+ 

AP 

Aran  of  Spill 

•f 

R. 

Solar  inaolation 

+ 

+ 

0 

Wind  speed 

+ 

L 

Spill  Depth 

+ 

ta 

Air  Temperature 

+ 

n 

Ground  Roughness 
Factor 

+ 

% 

Emineivity  of  the 
Atmosphere 

+ 

>KP  » 1, 

2,  3 • 1 - N2H4,  2 - MMH, 

3 - UDMH 

SECTION  IV 


ATMOSPHERIC  DISPERSION  MODEL 


I . Computational  rrocodure 

A elmple,  point  source  Gaussian  plume  model  was  selected  to  describe 
the  dispersion  of  propellant  vapors  released  into  the  atmosphere.  A com- 
parison of  Gaussian  models , Non-Gaussian  models,  including  the  tempera- 
ture-difference  model  used  for  safety  prediction  at  AF  missile  sites, 
box  models  and  the  gradient  transport  models  have  shown  the  Gaussian 
calculation  technique  to  be  the  method  of  choice  for  generalised 
usage.  Experimental  diffusion  studies  by  Cramer,  Haugen,  Paaquill  and 
others  have  confirmed  the  virtues  of  the  Gaussian  distribution  funation 
and  have  concluded  that  "tha  Gaussian  plume  formula  should  have  a wide 
area  of  practical  applicability  in  the  atmosphere"  (Reference  S) . 

The  Environmental  Protection  Agenoy  (EPA)  has  concluded  that 
"Gaussian  models  are  generally  state-of-the-art  techniques  for  estimat- 
ing the  impact  of  non-reactive  pollutants)  (Reference  6).  Furthermore, 
Gaussian  models  represent  the  best  choice  for  most  point  source  evalua- 
tions and  have  been  found  to  provide  reasonable  concentration  estimates 
(accuracy  by  a factor  of  2)  in  flat  or  gantly  rolling  terrain  (Refer- 
ence 6) . Although  Gaussian  models  have  been  praised  in  terms  of 
simplicity,  flexibility  and  good  correlations  between  predicted  and 
measured  downwind  concentration  values,  they  are  subject  to  the  follow- 
ing assumptions  and  limitations! 

a.  The  terrain  in  the  region  of  interest  is  either  flat  or 
gantly  rolling. 

b.  No  topographic  obstructions  occur  in  ths  vicinity  of 
the  source. 

c.  The  wind  speed  direction  and  other  meteorological  con- 
ditions ramain  uniform  and  oonatant  with  height. 

d.  Aerodynamic  downwash  does  not  occur. 


The  following  equation  is  used  to  calculate  ground  level  concen- 
trations for  a point  source,  continuous  release  emission  into  the 
atmosphere i 


C(xy) 


10° 

3600  ii  o o_  U 

y 1 


(3) 
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where 


t(*y> 

Cm 

0 


*>  concentration  at  coordinates  x,  y,  mg/m3 
= mass  transfer  rate  of  propellant  into  the  air,  kg/hr 
=•  average  surface  wind  speed,  m/seo 


c , na  » standard  deviation  of  the  concentration  distribution 

y along  the  y (crosswind)  and  t (vertical  axes). 

*•  » 

x,  y,  z - orthogonal  coordinate  system  in  which  the  origin  is 
at  ground  level  and  x is  in  the  downwind  direction. 


3600 


conversion  factor  for  time 


This  equation  is  applicable  for  ground  level  sources,  such  as  a spill, 
with  no  effective  plume  rise. 

The  constants  usad  to  calculate  the  aigma  y and  aigma  s 
coefficients  are  illustrated  in  Tables  3 and  4 for  various  stability 
categories  and  downwind  distanoss.  Thass  coefficients  are  currently  used 
in  various  SPA  models  (Reference  6)  and  ere  consistent  with  Turner's 
Workbook.  (Rsfsronca  6) . 


Equation  3 is  applicable  during  oonditione  whan  the  plume  is  un- 
restricted in  its  vertical  oxpansion.  During  inversions,  the  vertical 
expansion  of  the  plums  is  restricted  when  it  reaches  the  lower  bate  of 
a stable  air  mass.  This  type  of  inversion  occurs  most  frequently  during 
the  late  night  and  oarly  morning  hours.  The  vertical  concentration  pro- 
file becomes  uniform  and  is  reprssentsd  by  ths  following  equation  (Ref- 
erence 7) « 


C(xc) 


106  Qm 

3600  yj2n  oyHM  0 


exp- 


(4) 


where  H ■ height  from  ground  level  to  the  base  of  the  stable  layer  in 
meters  M 

A reasonable  assumption  is  that  the  vartioal  concentration  profile  begins 
to  be  affected  by  the  lid  at  that  diatanca  (XL)  downwind  from  tha  souroa 
whsra  concentration  at  ths  stable  layer  boundary  la  ena  tanth  of  tha 
plume  centerline  concentration.  XL  is  that  downwind  distanoa  from  tha 
source  where  cB  equals  1W2.15.  it  can  be  easily  determined  by  using 
figure  4 end  noting  the  downwind  distance  where  equals  Hm/2.15  (Raf- 
aranaa  7) . At  2Xr  and  graatar  ona  can  assuma  the "vertical  concentra- 
tion gradient  to  be  uniform  and  aquation  4 to  ba  applicable, 
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TABLE  3.  CONSTANTS  USED  TO  CALCULATE  Jy  IN  BPA  DEVELOPED 
DISPERSION  PROGRAMS  (REFERENCE  6) 


o a 1000  x. sin_ 0 » 465,12  x tan  0 
y 2.15  ooa  0 

x “ downwind  diatanoo  in  kilometers 

Stability 
(Categories*) 

A 

B 

C 

D 

0 

F 


Value  of  0 
(degress) 

24.167-2.5334  log#x 

18. 333- 1,8096  log#x 
12.5-1.0857  loggx 

8. 333- . 72382  log#x 
6. 25-. 54287  loggX 
4.1667-.36191  log,x 


♦stability  Categories  are  baaed  on  data  in  Turner's  Workbook  of 
Atmoaphorlo  Dispersion  Eatlmatea  (Reference  7) . 
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TABLE  4.  CONSTANTS  USED  TO  CALCULATE  a IN  EPA  DEVELOPED  PROGRAMS 

E 

(REPERENCE  6) 

CB  » ax*1 


Stability 
(Catagoriaa* } 

Downwind 

Diatanoo  (x)  (kilomator) 

Valuta  tor 

A 

o_  (m) 

* b 

A 

<0.1 

122. B 

0.9447 

0.1  - 0.15 

158.0 

1.0542 

0.15  - 0.2 

170.22 

1.0932 

0.2  > 0.25 

179.52 

1.1262 

0.25  - 0.3 

217.41 

1.2644 

0.3  - 0.4 

25B.B9 

1.4094 

0.4  - 0.5 

346.75 

1.7283 

0.5  - 3.11 

>3.11 

453.85 

** 

2.1166 

B 

<0.2 

90.673 

0.93198 

0.2  - 0.4 

98.483 

0.98332 

0.4  - 35.0 
>35 

109.30 

** 

1.0971 

c 

<123.0 

>123.0 

61.141 

• * 

0.91465 

D 

<U.  3 

34.459 

0.86974 

0.3  - 1.0 

32.093 

0.81066 

1.0  - 3.0 

32.093 

0.64403 

3.0  - 10.0 

33.504 

0.60486 

10,0  - 30.0 

36.650 

0.56589 

>30.0 

44.053 

0.51179 

E 

<0.1 

24.260 

0.83660 

0.1  - 0.3 

23.331 

0.81956 

0.3  - 1.0 

21.628 

0.75660 

1.0  - 2.0 

21.628 

0.63077 

2.0  - 4.0 

22.534 

0.57154 

4,0  - 10.0 

24.703 

0.50527 

10.0  - 20.0 

26.970 

0.46714 

20,0  - 40.0 

35.420 

0.37615 

>40 

47.618 

0.29592 

**Tha  maximum  calaulatad  valua  for  o is  5000  matars. 

♦stability  catagoriaa  art  baaad  on  dlta  in  Turnar'a  Workbook  of 
fttBQlPhirlc  Plaperalon  Eatlmataa  (Rofaronoo  7). 
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TABLE  4 (Continund) 


Stability 
(Catoqorif  *) 

V 


I 

i 


l 


Downwind  Valuaa  foe  o^  (m) 

Dlatancw  (x)  (kllomatara)  _a b 


<0.3 

15.209 

0.81538 

0.2  - 

0.7 

14.457 

0.78407 

0.7  - 

1.0 

13.953 

0.68463 

1.0  - 

2.0 

13.933 

0.63227 

2.0  - 

3.0 

14.B23 

0.34503 

3.0  - 

7.0 

16.187 

0.46490 

7.0  - 

15.0 

17.836 

0.41307 

13.0  - 

30,0 

:2.6S1 

0.32681 

30.0  - 

60.0 

27.074 

0.27436 

>60. 

0 

34.219 

0.21716 

I 

I 
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Figur*  4.  Vtrtiotl  Diip*r»ion  Coiffioiant  ab  a Function  of 
Downwind  DiBtAnoA  from  th«  Houma. 
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2.  Dispersion  Examples 


A sample  dispersion  problem  is  presented  to  demonstrate  the 
solution  technique  for  calculating  the  downwind,  ground  level,  center- 
line concentration  and  the  crosswind  distance  to  a specific  concentra- 
tion limit.  The  spill  sizes  and  specific  evaporation  rates  assumed 
were  taken  from  Figures  1-3,  for  a I5°c  ambient  air  temperature.  The 
ground  level  centerline  plume  concentration  is  calculated  using  the 
following  Gaussian  aquation i 

C , L°6  9M  (5) 

x ' *3600  "it""  7 U 

Equation  5 assumes  a high  mixing  layer  boundary  and  therefore,  minimum 
vertical  restriction  of  the  plume.  A 3m/sec  windapsed  and  Claes  B 
stability  oategory  representing  strong  to  moderate  solar  radiation, 
were  selected  for  this  example  (Reference  7) . Table  S lists  EPA 
calculated  horizontal  and  vertical  dispersion  coefficients  applicable 
to  this  problem  based  on  the  constants  described  in  Tables  3 and  4. 
Table  6 lists  the  center-line  concentrations  computed  for  a specific 
spill  size  using  Equation  5.  The  following  formula  is  used  to  calcu- 
late the  orosswind  distance  to  a specific  concentration  isoplethi 


where  y » crosswind  distance  from  the  plume  aenterline  to 
the  STPL 

Cxy  <*  crosswind  concentration  isopleth  ■ 7 mg/m3  for 

n2h4 

Cv  > downwind  oentsrline  concentration,  mg/m3 

! 

This  equation  ia  useful  in  determining  the  width  of  the  haserd  corridor 
associated  with  any  salaoted  limiting  concentration.  For  hydrasina  a 
reasonable  limiting  value  ia  7 mg/m3,  the  1 Hour  Short-Term  Public  Limit 
(stpl)  (Reference  8).  The  8TPL  for  MMH  and  UDMH  are  2,8  and  38.2  mg/m3 
respeotively.  Thee#  values  era  ussd  in  Figuras  5-7  and  17-19,  Table  7 
lists  tha  downwind  diatanoa,  oentsrline  ocnosntration,  and  orosswind 
distsnos  to  ths  7 mg/m3  isopleth.  Figures  5,  6 and  7 illustrate  the 
STPL  ieopleths  for  hydrazine,  MMH,  end  UDMH,  based  on  spill  conditions 
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TABLE  5 , 

Calculated  Standard  Deviations 

Used 

in  Sample  Dispersion  Problem* 

Downwind 

Horizontal  and  Vartic 

Distance 

Standard  Deviations 

(m) 

fa) 

if 

25 

5.42 

2.62 

50 

10.2 

5.18 

75 

14.8 

7.71 

100 

19.3 

10.6 

150 

27.9 

15.5 

200 

36.2 

20.2 

225 

40.2 

22.7 

250 

44.3 

25.2 

275 

48.3 

27.7 

300 

52.2 

30.1 

325 

56.1 

32.6 

350 

60.0 

35.1 

375 

63.9 

37.5 

400 

67.7 

40.0 

425 

71.5 

42.7 

450 

75.3 

45.5 

475 

79.0 

48.3 

500 

82.8 

51.1 

550 

90.2 

56.7 

600 

97.5 

62.4 

650 

104.8 

68.1 

700 

112.0 

73.9 

750 

119.1 

79.7 

800 

126.2 

85.6 

850 

133.3 

91.5 

900 

140.3 

97.4 

950 

147.2 

103. 

1000 

154.1 

109. 

1100 

167,8 

121.3 

1200 

181.4 

133.5 

1300 

194.8 

145.8 

1400 

208.1 

158.1 

1500 

221.3 

170.5 

1550 

227.9 

176.8 

1600 

234.4 

183.0 

♦Coefficients  calculated  using  equation!  in  Tables  10  and  11, 
Stability  Category  B 


TABLE  6.  Downwind,  Ground-Level  Centerline  Concentration 
for  Sample  Dispereion  Problem* 


Downwind  Distance 

Cx 

(m) 

mg/m' 

100 

' 

134 

150 

63 

200 

37 

225 

30 

250 

25 

275 

20 

300 

17 

325 

15 

350 

13 

375 

11 

400 

10 

425 

9 

450 

a 

475 

7 

500 

6 

*Aseumas  a 36000  liter  spill,  evaporating  at  930  kg/hr  (See  Figure  1. ) 
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TABLE 

7.  Croaawind  Distance  Valuea  To  Tha 

: 1 

7 mg/m3  Hazard  Corridor* 

1 

Downwind 

Downwind,  Ground- 

Croaawind 

Diatanca 

Level,  Center line 

Diatanoe  1 

(m) 

Concentration 

(mg/m3) 

y 

(m) 

100 

134 

47 

ISO 

63 

SB 

200 

37 

66 

225 

30 

69 

250 

25 

71 

275 

20 

70  1 

300 

17 

69.5 

325 

15 

69 

350 

13 

67  j 

375 

11 

61  !■ 

400 

10 

57 

425 

9 

51  j 

450 

B 

39 

475 

7 

0 

*Ha*ard  Corridor  defined  by  tha  Short-Term  Public  Limit  for  Hydraaina 
(Reference  7).  A 3600  litar  apill,  evaporating  at  930  kg/hr  is 
aaaumad.  Wind  Spaad  • 3 m/a,  Stability  B and  unraatriatad  mixing  daptha. 


CONDITIONS 


Figure  6.  Illustration  of  the  2.8  ag/m3  Short-Term  Public  Limit  Hazard  Cooridor  for  HMH. 


listed  in  Figures  1-3,  Ta  ■ 15°C,  end  unreetrioted  vertical  expansion 
of  the  plume.  Figures  8-16  present  the  hydrasine  oenter-line  ground- 
level  concentrations  for  200,  20,000  and  36,000  liter  spills  computed 
for  0,  15,  and  30°C  ambient  air  temperatures.  These  volumes  may  be 
associated  with  an  accidental  spill  for  a drum,  trailer  or  railcar, 
respectively.  Finally,  the  downwind  distance  to  the  BTPL,  as  a function 
of  evaporation  rata,  for  any  stability  category,  is  presented  in  Figures 
17-19.  These  graphs  represent  a quick  means  for  determining  the  down- 
wind distance  that  must  ba  evacuated  to  avoid  exposures  in  excess  of 
the  STPL.  These  curves  ware  calculated  based  on  a mean  wind  speed  of 
3 m/e  and  on  the  STPL  for  the  specific  propellant.  To  use  these  curves 
with  other  wind  speeds  or  limiting  concentrations,  the  evaporation  rata 
computed  for  a given  spill  situation  must  be  modified  as  follows i 

Q2  “ Qj  / STPL 

\cxy 


whore 

g2  ■ evaporation  rate  for  use  in  Figures  17,  16  and  19. 

Ql  - original  evaporation  estimate  from  program 

cxy  “ concentration  limit,  mg/m3 

U - mean  wind  apeed,  m/seo 


For  example,  aseume  a hydrasine  spill  for  the  following  conditions! 


(a) 

Ql  - 

1000  kg/hr 

(b) 

u - 

6 m/aeo 

(c) 

CXy 

2 mg/m3 

<d) 

Stability  Category 

Since  Figure  17  was  developed  for  an  STPL  concentration  of  7 mg/m3 
and  a wind  speed  of  3m/saci  the  amission  sstimets  of  1000  kg/hr 
must  be  modified  aa  follows! 

Q,  ■ 1000  MM  (8) 

w) 

g2  - 1750  kg/hr 


Q2  oen  be  directly  applied  in  Figure  17  to  estimate  e downwind  evacuation 
distance  of  2.2  km  to  the  2 mg/m3  concentration  limit. 
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IND  CONCENTRATION  - mg/m 
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IND  CONCENTRATION  - mg/m 


Figura  11,  Downwind,  Qround- 
a 200  Litar  MMH 
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CONDITIONS 

A Qm  * $00,1700, itQO  Af/hr 

n*  r4«o,/a,so*c 

I,  CLASS  • S TMILITY 
C,  0 •im/t 

0.  IWLL  OfATH  >t  tM 


.T4»  JO#C 


,t4*is*c 


-r4  ■ o*c 


DOWNWIND  DISTANCE  - m 

Figura  12.  Downwind,  Ground-Lava 1 Cantarlina  Concentration*  tor 
* 20,000  Li tar  MMH  Trailer  flpiil. 


c < 


CONDITIONS 


1 Figura  16.  Downwind , Ground-Lava 1 Cantarlina  Conoantrationa  for 
a 36,000  Litar  UDMH  Spill. 
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EVAPORATION  RATE  - kg/kr 


3.  Comparison  of  Hand  Calculated  Downwind  Dispersion  Data  with 
EPA  PTDIS  MOdsl 


Thu  downwind  concentration  data  graphically  illustrated  in  the 
previous  aeotion  ware  hand- calculated  using  dispersion  algorithms 
from  Turner's  Workbook  of  Atmospheric  Dispersion  Estimates  (Reference  7). 
This  hand-oaloulatad  data  was  then  compared  to  aomputer  generated 
concentration  data,  for  identical  spill  and  meteorological  conditions, 
using  the  SPA  point-distance  (PTDIS)  dispersion  program  (Reference  9) . 

The  efa  aomputer  program  verified  the  accuracy  of  the  hand- calculated 
dispersion  data  and  the  resulting  graphs. 

PTDIS  is  one  of  three  steady-state  Qaussian  plume  point  source 
models  that  have  recently  been  added  to  SPA'a  UNAMAP  (User's  Network 
for  Applied  Modeling  of  Air  Pollution)  system.  The  model  determines 
the  variation  of  ground-level  concentration  with  downwind  distanoe 
and  the  croaswind  distance  to  a selected  isopleth  concentration, 

Program  input  parameters  include  source  strength,  effective  height  of 
emission,  physical  steak  height,  stack  gas  temperature  stack  volume 
flow,  gas  velocity,  ambient  air  temperature,  stability  olass,  wind 
speed  and  mixing  height.  Concentration  values  for  up  to  50  downwind 
distances  may  be  computed  for  any  souraa  strength.  An  isopleth 
option  may  bo  aalled  which  calculates  the  half-width  to  a specified 
exposure  limit.  The  model  assumes  no  topographic  obstructions 
in  the  viainity  of  the  source  and  that  the  terrain  is  flat  or  gently 
rolling. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


Propellant  evaporation  and  dispersion  models  have  bean  developed 
for  hydrasine,  MMH  and  UDMH  propellant  ground  spills.  The  evaporation 
model  computes  the  rate  of  propellant  evaporation  as  a function  of 
ground  temperature!  solar  insolation,  air  temperature,  wind  speed 
and  spill  dimensions.  This  model  will  significantly  improve  pre- 
dictions of  hasard  sones  resulting  from  propellant  spills.  Results 
can  be  an  order  of  magnitude  different  than  those  obtained  using 
models  that  do  not  consider  heat  gain  and  losses  to  the  evaporating 
pool. 


Sensitivity  analyses  of  ths  evaporation  input  parameters  indicate 
that  spill  area,  ground  temperature,  ground  roughness  and  propellant 
type  are  important  factors  for  consideration.  Medium  sensitive  evapora- 
tion input  parameters  include  wind  speed,  air  temperature,  ground 
temperature  (when  Ta  is  less  than  T_)  and  mid-day  solar  insolation. 

Spill  depth  and  atmospheric  emiasivlty  were  found  to  be  the  least 
sensitive  in  the  evaporation  program. 

A simple  Qaussian  dispersion  model  is  presented  and  applied  to  a 
■ample  problem  to  calculate  the  peak  downwind,  ground- level  concentra- 
tion and  the  crosswind  distance  to  the  7,  2.8  and  38.2  mg/tir  8TPL.  This 
resulted  in  hasard  aorridors  reaching  480,  13B0,  and  680  meters  down- 
wind of  a 36,000  liter  railcar  spill  at  an  ambient  air  temperature  of 
1S°C  for  N]H4,  MMH  and  UDMH  respectively.  Downwind  hand-calculated  con- 
centration predictions  were  verified  and  found  to  agree  with  the  BPA 
point-distance  (PTDIS)  computer  model. 

Three  follow-on  projects  are  recommended.  The  evaporation  of 
propellant  mixtures  needs  to  be  addressed.  Titans  use  a 50  percent 
mixture  of  hydracine  and  UDMH  (Aarosane  50)  and  the  F-16  uses  a 70 
peroent  mixture  of  hydrasine  in  water.  Comparison  with  experimental 
laboratory  data  indicate  that  the  model  predicts  only  the  initial, 
first  hour  average  evaporation  rate.  The  parameters  that  describe 
the  transient  propellant  evaporation  rate  as  a function  of  time  and 
humidity,  need  to  be  defined.  Finally,  atmospherio  photochemical 
reaotion  or  decomposition  of  hydrasine  propellants  require  further 
investigation. 
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APPENDIX  A 


EMISSION  COMPUTATION 

Tho  avaporation  of  a liquid  pool  can  be  related  to  the  concentration 
driving  force  and  the  mass  tramfer  rate  constant  by  the  following 
equation i 

QM  “ km  <«1  * CA>  <A-1) 

where 

Qm  • mate  tranafer  rate,  kmoles/hr-m2 

k ■ vapor  phaae  maaa  tranafer  rate  coefficient,  m/hr 
m 

C ■ propellant  vapor  concentration  at  the  liquid  pool 
interface,  kmolea/m3 

m propellant  vapor  aonoentration  in  the  bulk  atmoaphere, 
aeaumed  aero. 

The  confutation  aaaumea  ideal  gas  behavior  of  the  propellant  vapor 
aoroea  the  diffusing  film.  The  concentration  term  (CJ.)  can  be  expreaaed, 
according  to  the  ideal  gaa  law,  in  vapor  preaaure  as  followai 

M/V  - Ca  - Pv/RTp  (A-2) 

and 

°M  " kmVR<Tp  + 273>  (A-D 

where 

Pv  * v*P°r  preaaure  of  the  apilled  liquid  kN/m2  ■ kPa 

R ■ universal  gaa  constant,  8,314  kPa-m3/kmole  • °K 

Tp  ■ temperature  of  the  evaporating  propellant  liquid,  °C 

M « number  of  moles  of  propellant,  kmolas 

V - volume  of  propellant  vapor,  m3 

The  mass  transfer  coefficient  (k^)  is  computed  from  the  following 
equationi  (Reference  3) 

X.  . CO  X (».«) 


4a 


where 


■ mass  tranafar  coefficient,  m/hr 

C - dimanaionlaaa  constant  and  a function  of  tha  Schmidt 
numbar 

U ■ wind  speed  at  a height  of  10m,  m/hr 
n « a function  of  tho  utmonpheric  lapse  rata 
X - diameter  of  spill,  m 

For  average  atmospheric  conditions,  a value  of  n ■ 0.25  is 
reasonable  (Reference  3) . The  constant  C is  calculated  for  any  pro- 
pellant using  ths  following  equation  (Reference  3)» 

C - 0.0392  So  “°*67  (A-5) 

where 

so  ■ Schmidt  number  - urM/pjrMDv 

Mpm  * viscosity  of  the  air-propellant  mixture  in  the  stagnant 
film,  g/am-seo 

P.M  - density  of  the  air-propellent  mixture  in  the  stagnant 
J’  film,  g/omr 

Dy  > diffusivity,  oma/sec 


The  sahmidt  number  is  e dimensionless  numbor  and  a function  of  the 
viscosity  and  density  of  the  air-hydrasine  mixture  in  the  stagnant  A'ilm 
and  the  diffusivity  of  the  propellant  fuel  in  air.  The  viscosity  of  the 
air-hydrasine  vapor  mixture  ia  based  on  the  pure  component  visocsities 
of  the  air  and  propellant  vapors  in  the  diffusing  film.  The  viscosity 
of  a pure  vapor  is  oalculatod  from  the  following  equation i (Reference  10) 


2V  (MWp)1/*  <Tf  + 273)  V2 

Vo  2/3  (Tr  + 1.47  Tb  + 674731)  10* 


(R*6) 


where 


uF))  - vieoouit.y  of  the  pure  qua  component,  g/cm-sec  or  poise 
MWp  ■ molecular  weight  of  the  propellant. 
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T » tamparatura  of  tha  stagnant  vapor  film 
- (TA  + Tp)/2,  oc 

Tft  ■ tamparatura  of  tha  air,  °c 

T^j  ■ boiling  point  tamparatura  of  propallant  liquid,  °C 

V0  ■ volume  of  tha  propallant  liquid  at  ita  boiling  point, 
cc/g-mole 

Tha  viscosity  of  a propallant  vapor  as  a function  of  tamparatura  can 
ba  raprasantad  by  tha  following  genoral  aquationi  (Raforsnca  10) 

«<v  2/1  <»-7> 

Tp  + 1.47  Tb 


whara 

K ■ constant  for  any  specific  gas 
Tb  « boiling  point  tamparatura,  °X 


A ganaral  aquation  for  tha  visooaity  of  air  as  a function  of  tampers- 
tura  oan  ba  computed  by  looking  up  a rafsranca  viscosity  value  at 
tamparatura  T-,  and  calculating  tha  constant  X from  aquation  A-7.  based 
on  this,  the  following  aquation  for  viscosity  of  air  was  developed i 


u 


FA 


1 .45  (10“5) 


T 3/2 
1F 

T_  4-  116 
F 


(A-8) 


The  viscosity  of  t)ia  vapor  mixture  in  the  stagnant  film  (upM)  is 
computed  from  the  following  aquation i (Reference  10) 


V Upp  (MWp)1/2  + (1  - V)  <mma)1/2 

y ■ 

Y(MW  )l/a  + (1  ..  Y)  (MWa)1/2 


(A-9) 
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where 


mpp,  yFA 

MW  , MW. 
P A 


viscosity  of  ths  sir-propellant  mixture,  poise 

avorage  mole  fraction  of  ths  propellant  vapor 
■ pv/2pt 

vapor  prasauru  of  the  pure  propellant  vapor,  kPa 

barometric  pressure,  ■ 101.3  kpa 

viscosity  of  the  pure  film  components,  propellant, 
air  respaotively,  poise 

molecular  weight  of  propellant,  air,  respectively 


To  solve  equation  A-9,  the  vapor  pressures  (Pv)  of  the  pure  gas 
components,  as  a function  of  the  pool  temperature  (Tp) , is  required. 
The  following  equations,  modified  for  unit  consistency,  are  pre- 
sented i (Reference  11) 


log  P..  - -7. .18113  - 653.880  + 0.047914  (T„  + 273)  (A-10) 

V™  P 

-4.90860  <10”5)  <Tp  + 273)a 


MMH 

log  Pv  - 6.23648  - 1104,571  - 152227.6  (A-U) 

T + 273  (T  + 373) 2 
P P 


UDMH 

log  Fy 


5.8*068 


875.89  - 140001.1 

tJ  '+  273  TtTT173)2 
p p 


(A-12) 


The  density  (ppM)  of  the  sir-propellant  mixture  in  the  stagnant  film 
must  alto  be  determined  to  evaluate  the  Schmidt  number.  The  density  of 
the  gas  mixture  is  represented  by  the  following  equation! 
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(A— 13) 


MW  Prr 


R'  T, 


where 

MW  - molecular  weight  of  the  propellant-air  mixture. 

- YMWp+<l-Y)MWA 

R'  - universal  gee  constant,  8314  kPa-em3/gmole  - °K 

Finally,  an  eetimate  of  the  diffusivity,  D , across  the  gas  film 
is  determined  from  the  following  equation i (Reference  12) 

°v  " °‘ 0010583  [V  ( 1/MWa  4 1/MWp)]  l/a  (A_14) 

Pm  *2 
T A-P  AP 

where 

Dv  - diffueivity,  cm2/sec 
rA-P  " <rA  + V'2 

r A,  rp  ■ collision  diameter  of  the  air-propellant  molecules, 
v Angstroms 

flA_p  ■ collision  integral  and  a function  of  kTj,/r.A_p 
k - Boltsman  constant,  1,38(10”6)  ergs/molecula-°K 
CA-F  “ energy  of  molecular  interaotion,  ergs 

The  collision  integral,  f)A_p,  oan  be  determined  by  calculating  first 
the  (cA_p/k)  force  constant  fur  each  gas  pair  as  followsi  (Reference  13) 


kT. 


A-P 


nv  2 


(A-13) 
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for  air,  i;  /k  “ lJV°K.  For  the  propallantB,  e / k was  calculated 
from  the  following  nquatiom  (Reference  13)  p 


Lp/k  * 1.15(Tb  + 273)  <A-1«) 

where  Tb  is  the  boiling  point  temperature  of  the  propellant  in  °c. 

Based  on  tabular  values  of  kTpA&_p  for  ths  temperature  range  of 
-31°c,  to  59°c  a general  regression  formula  for  tha  oolliaion  integral 
(fiA_p)  as  a function  of  temperature  was  developed) 


log  f!A_p  - -0.43  log 


/ k Tp\  + 0.15 

\1*-*  ) 


(A-17) 


The  following  specific  equations  were  developed  for  the  three 
propellants  under  consideration i 


N2H4 

log  f)A_p  ■ -0.43  log 

(Tp  + 273) 

+ 1.15 

(A- 18) 

MMH 

log  flA_p  * -0.43  log 

(Tp  + 273) 

+ 1.14 

(A-19) 

UDMH 

log  nA_p  - -0.43  log 

<Tr  + 273) 

+ 1.13 

(A-20) 

the  average  collision  diameter,  (rA  ) required  in 
collision  diameters  of  the  component  gas  species  ars 


To  calculate 
equation  A-14,  tha 
calculated  as  follows i (Reference  13) 


rp  - 1.18(V0)V3 


(A-21) 


where  VQ  ■ volume  of  the  liquid  at  the  normal  boiling  point, 
ac/g-mole 


The  collision  diameter  for  sir,  r«,  is  given  ee  3.617A  (Reference  13). 
v0  ia  calculated  by  dividing  the  density  of  the  propellant,  at  its  boiling 
point,  into  tha  molecular  weight  of  the  propellant,  The  density  of  a 
liquid  at  its  boiling  point  is  computed  by  the  following  equation) 
(Reference  10) 


& 


PLS  " PLl  / TC  " r 

\TC  - 


tli/ 


(A-22) 
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where 


p-a  - density  of  the  liquid  at  ite  boiling  point,  g/cc 
N2H4  - 0.911,  MMH  - 0.801,  UDMH  - 0.740 

PLI  - denaity  of  the  liquid  at  Ttj,  g/ee 

Tjj  ■ temperature  of  a liquid  at  its  boiling  point,  °C 

Tq  - critical  temperature  of  the  propellant,  °c 

Tl1  «*  any  propellant  temperature  where  denaity  is  known. 


Equation  A- 3 can  now  ba  solved,  assuming  a tentative  pool  tempera- 
ture, T_ . Tha  pool  temperature  of  an  evaporating  pool  may  not  equal  the 
temperature  of  the  air.  To  determine  the  steady-state  pool  temperature 
(Tp) , a heat  balance  analysis  of  tha  evaporating  pool  system  must  be 
performed.  The  following  parameters  define  the  heat-balance  equation) 
(Reference  3) 

Qg  + Qh  + Qs  + CA  - + He  (A-23 ( 


where 

Q ■ convective  heat  transfer  from  tha  noil  or  surfaaa  to  the 
pool,  j/hr 

Qh  " convective  heat  transfer  from  the  atmosphere,  J/hr 
Qs  * solar  insolation,  J/hr 

Qa  ■ radiative  heat  gain  from  the  atmosphere,  J/hr 
* radiative  heat  loss  from  the  liquid  pool,  J//hr 
Hg  * evaporative  heat  loas  from  the  pool,  J/hr 


Heat  transfer  from  the  ground  to  the  liquid  pool  (Qg)  is  based  on 
experimental  heat  transfer  studies  through  a horizontal  layer  bounded 
on  the  top  by  a cold  surfaaa  and  on  the  bottom  by  a heated  eurfaoe. 
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The  Rayleigh  number  <N  ) , (Graehof  number  (NQR)  time*  the  Prandtl 
Numbor  <Nt,R) ) is  correlated  with  the  Nussolt  number  for  five  distinct 
modes  of  hunt  transfer.  Tho  five  modern  arc  defined  as  follows t 


a.  100  < Njy^  < 1700 


N 


NU. 


- 1.  Heat  transfer  is  by  conduction 


only  and  the  fluid  remains  immobile. 


b.  1700  < Nv»  < 3000  NNu  ■ 0.0012  N^0,90  This  mode  defines 
the  oritical  condition  where  creeping  convection  begins.  The  fluid 
begins  to  circulate  and  heat  transfer  is  by  conduction  and  convaotion. 


a.  3000  < Nra  < 8000  NPR0-2  - 0,24  N,^0*25  This  heat 
transfer  mode  ie  defined  as  laminar  convaotion  and  is  described  by 
uniform  heat  flow  contours. 


d.  8000  NpR0*2  < Njy^  < .18000  NpR0,2  % - 0.3  Ncll°*16NpR0*21 

This  region  is  defined  as  tha  transitional  phase  where  heat  oonveotion 
change*  from  laminar  to  turbulent  convection. 


e.  Njy^  > 18000  NpR  0,2  Nhu  - 0.10NGR°'31NpR0,36 

This  region  is  aharaatsriasd  by  turbulsnt  convective  heat  transfer  from 
the  ground  to  the  liquid  pool. 


The  heat  tranafar  from  the  ground  to  the  liquid  pool  is  described 
by  the  following  equation i 


Qa  - hg  <Tg  - Tp)  Ap  (A-24) 

where 

h ■ heat  transfer  coefficient  from  ground  to  liquid 
’ pool,  J/m2,hr,0X 

TG,Tp  - tomperaturo  of  tho  ground  and  pool,  respectively, 

Ap  « area  of  tho  liquid  evaporating  pool,  m2 

S 

I 


1 


The  heat  transfer  coefficient  can  be  calculated  from  the  follow- 
ing aquatiom  (Rafaranca  14) 

v _ kL  A N0R  B NpR  C (a-25) 

9 L 


i 


i 

ii 


where 


1 \f 

i ' t 


L 

A,  B,  C 


N, 


’OR 


8 

L 

9 

Pl 

N, 


PR 


thermal  conductivity  of  tha  liquid  J/m*ur,{V 

ccinatanta  whoaa  value  dapand  upon  tha  fiva  distinct  modes 
of  haat  transfer  that  can  occur. 

12  2 

GJr**hof  numbtr  ■ 9 0 L p ¥/P  , 

h L 

ratio  of  volume  change  to  mean  volume  par  degree 
change  in  temperature,  1/°K 
liquid  layer  thickness,  maters 

gravitational  constant  « 1.27  x 108  m/hr2 

liquid  density,  kg/m3 

liquid  viscosity,  kg/m*hr 

Prandtl  number  ■ 


Depending  on  the  heat  transfer  mode,  values  of  A,  B,  C are  assigned. 
For  instance,  when  TG  < T_,  haat  transfer  to  the  ground  is  by  con- 
duction and  A » 1,  and  B > c • 0,  Equation  A-25  then  simplifies  as 
follows i 


(A-26) 


For  thie  case  the  heat  transfer  rate  is  inversely  proportional  to  the 
pool  depth.  When  Ta  > T and  heat  transfer  is  by  turbulent  convection, 
A * 0.1,  B - 0.31  and  C - 0.36.  For  this  case,  Equation  A-25  looks  as 
follows! 


k OlN  ® n 0.36 
kL  0,i  NQR  NPR 


(A-27) 
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l 


Slnoo  N„_  is  a function  of  l3,  pool  dapth  oancals  out  In  tha  above 
equation  and  becomes  an  lnaanaitive  parameter  in  the  heat  transfer 
equation. 

The  heat  transfer  from  the  atmosphere  to  the  evaporating  pool  is 
described  by  the  following  equation!  (Reference  3} 

Qh  - h(TA  - Tp)Ap  (A-28) 

h » heat  transfer  coefficient,  J/m2»hr,0K 

h is  further  defined  asi  (Reference  3) 

h " ViTM  CprM  iso/PR>0,67 


where 

dpm  « density  of  the  film  vapor  mixture,  g/m3.  For  unit 
cons latency,  multiply  the  density  in  g/co  by  10*. 

Crfm  " h,4t  capacity  of  film  vapor,  J/g,0K 

Equation  A-29  is  modified  as  shown  below  for  dimensional  oonslstenoyi 
h - 106  ppM  CpyM  <8c/PR)°*67  (A-30) 

The  heat  capacity  of  the  air-propellant  vapor  mixtura  in  tha  film 
ie  calculated  first  for  sach  purs  componsnt  vapor  as  follows*  (Refer- 


unce  13) 

w2«4 

(A-31) 

CPPP  " 

[^0.357  + 

7.919<10"4) (Ty  + 273) 

+ 2. 

44(10"7) (Ty  + 273)3]4.19 

MMH 

(A-32) 

CPPP  " 

8,49(10“2)  + 9.88(10"4) (Ty  + 

273) 

- 3.22 (10“7) (Ty  + 373) 2 J 4 

UDMH 

(A-33) 

CPPP  “ 

6.32(10' 

‘2)  + 1.087 (10-3) (Ty  + 

273) 

- 3.63 (10“7) <Ty  + 273) 2j 
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where 

CPFP  - haat  capacity  of  tha  pura  propellant  vapor,  J/g,0K 
Ty  ■ tamparatura  of  tha  film,  °C 


Air 


(A-34) 


CPFA  " [°*232  + 1.622(10“5)  (Tp  + 273)  + 3.96(10“®)  (Tp  + 273)2]  4.19 


where 

CPFA  « haat  capacity  of  air'  J/,g*°K 

Tha  haat  capaaity  of  tha  propellant-air  mixtura  ia  rapraaantad 
by  tha  following  aquation* 

r - V<CPFP)  WWp)  + (1  “*>  (CpyA)  (MWA)  ( 

(y)  +"(r-"y)"WA5 


Tha  Prandtl  numbar,  alao  raquirad  in  aquation  A-30,  ia  dafinad  by 
tha  following  expreaaioni 

P - CPFM  <W  3-6 <A“36> 

R C 

KFM 


whara 


CPFM 

^FM 

^FM 

3.6(10®) 


■ Prantl  number,  dimanaionlaaa 

■ haat  capacity  of  tha  vapor  film  mixtura,  j/g*°K 

■ viacoaity  of  tha  vapor  film  mixtura,  g/cm*aac 

« tharmal  conductivity  of  tha  film  mixtura,  J/hr *m,0K 
• csomrcraion  factor  to  changa  g/om*aac  to  g/m'hr 


Tha  tharmal  conductivity  (hyp)  of  tha  pura  propellant  vapor  ia 
confuted  from  tha  following  aquation*  (Rafarence  13) 
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where 


k * tnermal  conductivity  of  the  pure  propellant  vapor, 

J/hr*m*°K 

urp  * viaooaity  of  tha  propellant,  g/om*aec 

Cppp  ■ heat  oapaoity  of  vapor,  j/g*°K 

For  air,  tha  thermal  conductivity  (kp.)  ia  axpreaaed  by  an  aquation 
that  aaaumea  a atrnight-lina  relationahip  between  kpA  at  0°C  and  100°Ci 

kyA  - 87.104  + (Tp)  0.2679  (A-38) 


where 


kpA  ■ thermal  aonduotivity  of  air,  J/hr'm,0X 
Tp  ■ temperature  of  the  film,  °C 

Tha  thermal  oonduotivity  of  the  air-propellant  mixture  ia  rapraaentad 
by  tha  following  equation i (Reference  13) 


^FM 


Y kpp  (MWp) 1/3  + (1  - Y)  kpA  (MWa>  1/>3 

y (MM  ) 1/3  + (i.y)  (MW.)1/3 
P “ 


(A-39) 


The  heat  input  to  the  pool  due  to  aolar  radiation  ia  deaoribad 
by  the  following  equation!  (Reference  3) 


Qa  - (1-a)  R#Ap  (A-40) 

where 

Qa  ■ incident  aolar  radiation,  J/hr 

a « aurfaca  albedo,  the  radiation  that  ia  reflected  from  the  pool 
Ra  - aolar  radiation,  J/m2<hr 


The  albedo, is  tentatively  eeleated  ae  0,14  which  is  in  reasonable 
agreement  with  the  generally  accepted  value*  for  water  surfaces,  (Ref- 
erence 3)  . 

Thu  radiative  heat  input  to  t.ho  propellant  pool  due  to  oounter- 
radiation  from  the  atmosphere  ia  represented  by  the  following  equation! 
(Reference  3) 


Q.  - e S (Ta  + 373) 4A  (A-41) 

A * r 

where 

qa  • radiative  heat  gain  from  the  atmosphere,  J/hr 
e#  " emissivity  of  the  atmosphere 

s ■ Steffan-Boltsmann  constant,  2.042 (10“4)  j/m2 • hr • °K4 


The  emissivity  of  the  atmosphere  is  primarily  a function  of  the 
water  vapor  pressure  and  tentatively  assumed  constant  at  0.75  (Ref- 
erence 3) . 

The  radiative  heat  loss  from  the  pool  to  the  atmosphere  is  presented 
by  the  following  equation!  (Reference  3) 

Qb  ■ ep8(Tp  + 273)  4Ap  (A-42) 


where 


g_  • heat  loss  from  the  pool,  J/hr 

a 

e ■ emissivity  of  the  pool 


The  emissivity  for  water,  equal  to  0.95,  is  tentatively  used  until 
specific  values  for  the  hydrasina  family  can  be  obtained. 

finally,  the  heat  loss  from  the  pool  due  to  evaporation  oan  be 
expressed  by  ths  following  equation!  (Rsferenoo  3} 


H 


E 


K Ap  MW  Xp  PV 


(A-43) 
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where 


HB  - evaporative  heat  l0„  from  the  pool|  J/hr 

k ■ vapor  mass  transfer  rate  coefficient,  m/hr. 
m See  equation  A-4. 

X.  ■ enthalpy  of  vaporisation  at  the  equilibrium  pool 
temperature,  J/kg 


Using  Watson's  Correlation,  X can  be  computed  from  the  following 
expression!  (Reference  13) 


where 

Tc  ■ critioal  temperature  of  the  propellant,  °C 
X • reference  enthalpy  of  vaporisation  at  Tj_ 


The  critical  temperatures  for  the  three  liquid  propellants  are  as 
follows i N2H4  - 380°Ci  MMH  - 312°C;  UDMK  - 230°C.  The  reference 
enthalpy  of  vaporisation  (Xj.)  used  in  the  sample  calculations  era  as 
foliowsi  NaH4  - 1.36(10*)  J/kg  at  T,  - 230C|  MMH  - 0.7S(1O5)  J/kg 
at  T],  « 25°C)  UDMH  - 5.44<105)  J/kg  at  Tj  « 62.5°C. 


A value  of  Tp  is  assumed  and  equations  A-24  to  A-44  are  solved  to 
determine  if  the  equality  condition  set  by  equation  A-23  has  been  met, 
Successive  iterations  are  made  until  Tp  + 0.1  are  oaloulated.  This 
value  is  subsequently  used  in  equetionA-45  to  caloulate  the  mass 
flux  of  propellant  evaporating  into  the  atmosphere. 


QM  - kJnPv  (MW)  Ap/R  Tp 


(A— 45) 


where 

qh  ■ Emission  rate  from  the  pool,  kg/hr 
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APPENDIX  B 


EVAPORATION  PROGRAM  FLOW-CHART 


( START") 

X 


READ  DATA 
TAC.U,DSPILL1ASPILLi 
GN,RS,KPJG 
DP 


TP»TA 


COMPUTE  PROPELLANT  VAPOR 
AND  LIQUID  PROPERTIES  i 
COMPUTE  MASS  TRANSFER 


immmmcE 

COMPUTE  • TPT 
f(TPT)*0 


TPNEW 
s 0.001 


NO 


INCREMENT 

SELECTED 

VARIABLE 
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:TP  + FMC(TFT-TPi 


HYDRAZINE  EVAPORATION  PROGRAM 
TAC,U,DSPILL,ASPILL,GN,RS,KP,TG>OP 
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LIST  OF  SYMBOLS 


j t 

» ■ 


i 


COMPUTED  TERMS  IN  EVAP  PROGRAM 

Computed  Variable! Variable  Description 

Technical  Computer 
Report  Program 


A $ B | C f 

A,  B,  C, 

Conatanta  uaed  in  generalicad  heat  tranafer 
correlation  equation,  N^  ■ A NaRB  NpRc . 

B 

Beta 

The  ratio  of  volume  change  to  mean  volume  per 
degree  change  in  temperature,  1/°K 

Qs+Qa“he 

CONST 

Solar  inaolation  + radiative  heat  gain  from 
the  atmoaphere  - evaporative  heat  loea,  J/hr 

CL 

CPL 

Liquid  heat  capacity,  j/kg*°K 

CPFM 

CPFLM 

Heat  capacity  of  the  propellent-air  mixture 
j/g*°K 

CPFP 

CPP 

Heat  capacity  of  pure  propellant  vapor, 

J/g*°K 

PFM 

DENFM 

Deneity  of  the  air-propellant  mixture  -g/om3 

PL 

DENL 

Liquid  density  kg/m3 

DV 

DV 

Diffuaivity  of  propellant  gaa  aoroaa  the  gaa 
film,  cma/seo  - needed  for  Schmidt  No. 

ngr 

GR 

2 2 

Qraahof  Number  - gBAt  L3  p ^ / P L 

h 

H 

Heat  transfer  rate  coefficient,  J/ma*hr,0K 

CPFA 

HCA 

Heat  capacity  of  air,  J/g,wK 

\ 

HG 

Heat  transfer  rate  coefficient,  from  the  ground 
to  the  liquid  pool,  J/ma,hr,0K 

UA-P 

OMEGA 

Collision  integral  for  diffusion,  function  of 
kTAA_p 

PR 

PR 

Prandt's  number  for  the  vapor  phase  ■ 

CFFM  uFM'/kPM 

npr 

PRL 

Prandtl  No.  for  the  liquid  phase  ■ 

Nra 

RA 

Rayleigh  No.  ■>  nqr  * NPR 

"Visa—— Tp 
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LIST  OF  SYMBOLS  (continued) 


Computed  Variables Variable  Daaorlgtlgn 


Technical  Computer 
Report  Program 


rA_P 

RBAR 

Average  collieion  diameter  of  the  air  pro- 
pellant moleoulee,  A° 

sc 

SC 

Schmidt  number,  dimenaionleeB  - u/oDy 

ta 

TA 

Temperature  of  Air,  °C 

tf 

TF 

Average  temperature  vapor  film,  °C  - 
<Tp  + TA)/a 

kFM 

TCOFM 

Thermal  conductivity  of  the  propel lent- air 
mixture 

kFA 

TCONA 

Thermal  conductivity  of  air,  J/hr*m*°K 

kL 

TCONL 

Liquid  thermal  conductivity,  J/m*hr*°K 

krp 

TCONP 

Thermal  conductivity  of  the  pure  propellant 
vapor  j/hr*m*°K 

,pl' 

TP 

Computed  temperature  of  the  pool,  °K 

- 

TFT 

Temporary  pool  temperature  calculated  by 
ueing  Raphoon-Newton  method 

Qh  + 

Qg 

TPTi 

Convective  heat  transfers  from  (to)  the 
atmosphere  + from  the  soil,  J/hr 

qe  + 

+ 

TPT 

2 

Radiative  heat  loss  from  the  liquid  pool  + 

Qa’ 

H 

a 

CONST  J/hr 

IR 

TRX 

Mass  transfer  rate  coefficient,  m/hr 

VAPRT 

Emission  rate  from  the  liquid  evaporating 
pool,  kg/hr 

11  FA 

VISA 

Viscosity  of  air,  poiee  » g/cm-sec 

mfm 

VISFM 

Viscosity  of  the  air-propellant  vapor 
mixture,  poise 

UL 

VISL 

Liquid  visooaity  Kg/m* hr 

,JFP 

VISP 

Viscosity  of  the  pure  gas,  component,  poise 
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LIST  OF  SYMBOLS  (continued) 


Computed  Variables 
Technical  Computer 
Report  Program 

Vc  VO 

Pv  VPP 

Y Y 


r 


P 


N 


NU 


Variable  Daeoriptlon 


volume  of  liquid  propellant  at  its  normal 
boiling  point,  oo/g-mole 


Vapor  pressure  of  the  pure  gaa  component*, 
kp*  ■ kN/m2 


Average  mole  fraction  of  the  pure  vapor 
component 


Collision  diameter  of  propellant  molecules  ■ 
1.18  v,,1'3  Angstroms 

Nuaaelt  Number  - hg  L/k^ 
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LIST  OF  SYMBOLS 
INPUT  DATA  TO  EVAP  PROGRAM 


Input  Symbol 


Technical 

Report 


Computer 

Report 


Variable  Desorigtlon 

Area  of  Spill 

Depth  of  pool  in  metera 

Diameter  of  spill,  m 

Ground  roughneee  factor  and 
temperature  profila  in  the 
atmosphere 

Liquid  propellant  type  to  be 
evaluated i KF  - 1,  2,  3 
1 - N2H4  3 - MMH , 3 - UDMH 

Solar  insolation  rate,  J/m2*hr 

Temperature  of  air,  °c 

Ground  temperature  at  10  am 
depth 

Wind  speed,  m/sec 


APPENDIX  C 


SENSITIVITY  ANALYSES 

Tablaa  c-1  through  C-ll  dapict  input  pirunttar  variation  and  tha 
resulting  avaporation  rataa  for  tha  following  initial  hydraiine  apill 
conditional 


a. 

Air  temperature 

15°C 

b. 

Wind  speed 

3 m/nec 

a. 

Spill  volume 

20  m3 

d. 

Spill  diameter 

31.66  m 

a. 

Spill  area 

787.40  ma 

f. 

around  roughness  factor 

0.25 

<3> 

Solar  rata 

3.1  (106)  J/ma-hr 

9. 

Solar  rata 

3.K106)  J/ma*hr 

h. 

Liquid  type 

hydraaina 

i. 

around  temperature  - 

air  tamparatura 

j. 

Depth  of  pool 

2,54  om 

Tabla  C-1  illuatrataa  tha  variation  in  evaporation  rataa  batwaan 
hydraaina,  MMH  and  UDMH.  In  general,  tha  avaporation  rataa  of 
udmh  and  MMH  ara  11  and  3 times  that  of  hydraaina.  Tabla  c-ll  shows 
that  tha  avaporativs  haat  loss  la  tha  major  factor  raduoing  tha  pool 
tamparatura,  Since  tha  avaporation  rata  ia  a function  of  tha  vapor 
praasura  of  tha  propallant  and  tha  mass  tranafar  coefficient,  tha  maao 
transfer  coefficient  must  bo  the  factor  increasing  tha  avaporation 
rata  while  raduoing  the  pool  tamparatura. 

Tabla  c-3  indicates  that  ground  temperature  ia  both  a high  and 
medium  aanaitivity  input  parameter  to  tha  evaporation  program, 
when  tha  ground  temperature  ia  below  the  equilibrium  pool  tamparatura, 
heat  tranafar  oaoura  only  by  aonduotion.  Table  C-3  indicates  that  a 
34  paroant  error  in  tha  avaporation  rata  occurs  for  every  10°C  error 
in  tha  ground  temperature  for  the  above  condition.  Whan  ground  tem- 
perature exoeeda  tha  pool  temperature,  heat  tranafar  usually  occurs 
by  turbulent  convection.  Tabla  C-3  indicates  that  a 55  paroant  error 
can  occur  for  this  condition.  Table  c-ll  shows  that  haat  tranafar  by 
aonvection  in  on  the  order  of  300  percent  greater  than  by  conduction. 
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Evaporation  rate  ia  directly  proportional  to  apill  area  and 
increases  non-linearly  with  the  spill  diameter.  „.Table  C-3  illus- 
trates the  variation  in  evaporation  rate  with  constant  increases 
in  spill  diameter.  Assuming  a 30  meter  (diameter)  spill,  a 33 
peroent  variation  in  the  spill  diameter  results  in  a 75  percent 
oversatimats  or  a 34  percent  underestimate  in  the  evaporation  rate. 

This  example  assumes  a 20  m3  trailer  spill. 

Table  04  indicates  that  an  error  in  the  insolation  rate  can 
produao  a proportional  error  in  the  evaporation  rate  for  high 
insolation  values.  Whan  the  inaolating  rate  is  low  (compared  to 
the  other  terms  in  the  heat  balanee  equation) , the  resultant  pool 
temperature  and  therefore  the  evaporation  rate  is  little  influenced. 

Table  C-4  shows  that  for  solar  insolation  below  0.8  MJ/nr>hr  (late 
afternoon  or  early  morning) , heating  effect  of  the  aun  is  insigni- 
cant.  When  insolation  is  above  0.8  MJ/m3,hr  (analysis  valid  only 
for  the  input  parameters  described  above) , the  heat  transfer  to  the 
ground  or  to  the  air  are  not  eufficient  to  balance  the  solar  heat 
input,  resulting  in  a higher  pool  temperature  and  increased  evapora- 
tion rate. 

Wind  speed  is  a medium  sensitive  parameter  in  the  propellant 
evaporation  program.  Table  05  shown  that  as  wind  speed  increases 
from  one  to  five  meters  per  second,  the  evaporation  rate  doubles. 

Table  Oil  indicates  that  the  primary  parameter  responsible  for  the 
increased  evaporation  rate  is  the  masa  transfer  coefficient  (km) , a 
function  of  velocity,  which  doubles  the  evaporative  heat  loss  from 
the  pool. 

Pool  depth  has  a negligible  effect  on  the  evaporation  rate. 

During  heat  transfer  by  conduction,  when  the  ground  temperature  ia  oooler 
than  the  equilibrium  pool  temperature,  a doubling  of  the  pool  depth 
resulted  in  a 12  percent  average  increase  in  the  evaporation  rate,  The 
increase  did  not  exceed  20  percent  for  any  specific  pool  depth.  During 
the  conductive  heat  transfer  phase,  the  heat  transfer  coefficient  to 
the  ground  reduces  to  ths  simple  expression) 


As  the  pool  depth  (&)  increases,  tu  decreases  and  the  pool  warms  up. 
During  the  convective  heat  transfer  phase,  h„  in  a more  complicated 
functioti  where  L is  also  in  the  numerator.  Thus,  when  pool  tempera- 
ture is  less  than  the  ground  temperature,  depth  is  not  a sensitive 
parameter  in  controlling  the  evaporation  rata.  Tabla  C-6 
illustrates  the  increase  in  pool  temperature  and  evaporation  rata  when 
pool  temperature  is  greater  than  the  ground  temperature. 
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The  ground  roughness  factor  (n) , a highly  sensitive  parameter, 
was  assumed  to  equal  0.25  for  average  atmospheric  condition 
(Reference  3) . Variation  of  the  pool  evaporation  rate  with  n is 
shown  in  Table  C-7.  In  general,  a 100  percent  change  in  n results 
in  a 43  percent  or  greater  change  in  the  evaporation  rate.  As  n 
increases,  the  velocity  term  (equation  A-4)  whioh  is  the  controlling 
parameter  that  affects  the  evaporation  rate,  decreases.  Although  Tp 
increases,  which  increases  the  heat  transfer  driving  force  between 
the  air  and  ground,  heat  loas  to  the  air  ia  limited  by  the  reduced 
velocity  and  heat  loss  to  the  ground  is  a function  of  the  pool  depth. 

Moat  attampts  to  datermina  the  range  of  n have  found  that 
0<n<l.  Ths  greatast  range  of  observed  valuta  ia  0.04<n<0.93  with 
n «•  0.23  during  conditions  of  adiabatic  lapse  rats  (Referenda  15). 

This  value  is  vary  close  to  ths  l/7th  power  law  familiar  in  turbulent 
pipe  flow  (Reference  15),  Table  C-8  presents  values  of  n for  various 
temperature  difference  data  in  the  etmoaphere  between  5 and  400  feet 
abovo  the  ground  (Reference  15). 

The  emissivity  of  the  atmosphere  (eA)  is  a function  of  ths  water 
vapor  pressure  in  the  atmosphere  end  assumed  es  0.75  in  this  report. 
Table  C-9  shows  that  ea  ia  not  important  in  controlling  ths  evapora- 
tion rats.  A reduction  in  ea  of  0,25  (33  percent)  results  in  a 10 
percent  reduction  in  ths  evaporation  rata. 

Variation  in  sir  temperature  from  0 to  40°c  results  in  a 100 
percent  increase  in  the  evaporation  rata,  in  general,  a 10°c  error 
in  ths  sir  temperature  results  in  s 20  percent  error  in  the  evaporation 
rats.  This  variable  ia  designated  as  a medium  aansitivs  parameter , 
Table  C-10  show*  the  variation  of  evaporation  rata  with  incraasing  air 
temperature . 
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TABLM  0-1.  VARtATTON  OP  EVAPORATION  RATE  WITH  PROPELLANT 


propoHum 

Pool  Temperature 

Evaporation  Rata 

°K 

kg/hr 

Hydraaina 

299.2 

538.5 

MMH 

293.2 

1704 

UDMH 

237.7 

5807 

TABLE  2-C.  VARIATION  OP  EVAPORATION  RATE  WITH  GROUND  TEMPERATURE 


Ground  Temperature 

Pool  Tamparatura 

Evaporation  Rata 

°K 

°K 

Rg/h* 

273 

293.3 

385.7 

283 

297.3 

483.8 

293 

301.1 

597.8 

303 

304.7 

727.0 

313 

312,  B 

<• 

1111.5 

316 

317.7 

1426.4 

TAB  Eli  03. 


VARIATION  OF  EVAPORATION  WITH  SPIU,  SIZE 


pJU  Diameter  Aria  Pool  Temperature  Evaporation  Rate 


m 

m- 

10 

7B.5 

20 

.1 14.2 

.10 

706.') 

40 

1256.6 

50 

1963.45 

°K 

kg/hr 

298.1 

57.8 

298.  0 

221 

299.  2 

485 

299,5 

848 

299.7 

1306 

TABLE  C-4.  VARIATION  OF 


Solar  Insolation 
MJ/m2hr 
0.1 
0.2 
0.4 
0.8 
1.6 
3.2 
6.4 


RATE  WITH  SOLAR  INSOLATION 


Pool  Tonjparatura 
°K 
288 
288 
288 
288 
292 
300 
314 


Evaporation  Rata 
kg/hr 
283 
283 
283 
283 
358 
569 
11B6 


TABLE  C-5.  VARIATION  OF  EVAPORATION  RATE  WITH  WIND  SPEED 


Wind  Spaed 

Pool  Tamparature 

Evaporation 

m/aac 

°K 

Rate  kg/hr 

1 

305 

319 

2 

302 

448 

3 

299 

539 

4 

298 

612 

5 

296 

674 

76 


TABLE  C-6.  VARIATION  OF  EVAPORATION  RATE  WITH  POOL  DEPTH 


Pool  Depth 

Pool  Temperature 

Evaporation  Rate 

Centimeters 

°K 

kg/hr 

0.15875 

290 

309 

0.3175 

291 

335 

0.635 

293 

381 

1.27 

296 

451 

2.54 

299 

538 

5.08 

302 

622 

10.16 

304 

683 

20.32 

305 

721 

TABLE  C-7.  VARIATION  OF  EVAPORATION  RATE  WITH  THE 
GROUND  ROUGHNESS  FACTOR 


Ground  Roughness  Pool  Temperature  Evaporation  Rate 

Factor  (n)  °K  kg/hr 


0.1 

269.1 

1220 

0.2 

295.8 

693 

0.3 

302.4 

419 

0.4 

307.4 

247 

0.5 

310.6 

140 

0.6 

312.4 

78 

0.7 

313.4 

44 

0.8 

313.9 

25 

0.9 

314.2 

15 

1.0 

314.2 

9 
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TABLE  C-9.  VARIATION  OF  EVAPORATION  RATE  WITH 
ATMOSPHERIC  EMISSIVITY 


Emissivity 

Pool  Temperature 

Evaporation  Rate 
kg/hr 

0.75 

299.2 

539 

0.65 

29B.S 

516 

0.55 

297.7 

494 

0.45 

296.6 

473 

0.  35 

296.1 

452 

f 


TABLE  C-10. 

VARIATION  OP  EVAPORATION  RATE  WITH 
AMBIENT  AIR  TEMPERATURE 

Air 

Pool 

Evaporation 

Temperature 

Temperature 

Rata 

°K 

°K 

kSf/hr 

273 

294.3 

408 

283 

297.7 

493 

293 

300.8 

588 

303 

303.9 

695 

313 

306,9 

818 

J 


SISATWre 


INITIAL  DISTRIBUTION 


HQ  USAF/LEEVP 

2 

Dat  1 ADTC/TST  (Taoh  Library) 

1 

HQ  USAF/SAFOI 

1 

Dat  1 ADTC/ECC 

2 

HQ  USAF/SOPA 

1 

SAMTEC/WE 

1 

HQ  ADCOM/SOPAP 

1 

DARCOM,  Field  Safety  Activity 

1 

NCEL,  Coda  25111 

1 

Southeast  Bnvmtl  Raoh  Lab 

1 

HQ  TAC/SQPB 

3 

DD-MED-41 

1 

HQ  USAFE/Surgeon 

1 

AFOSR/NC 

1 

HQ  USAFE/DEPV 

1 

AWS/LNP 

1 

HQ  AFXSC/SQMS 

1 

HQ  SAC/DOW 

} 

HQ  USAFA/DEV 

1 

AFIUPL/WE 

\ 

Toxic  Material*  Info  Ctr 

1 

Dat  1 ADTC/r\Y 

a 

AFOSR/Lif • Soiancaa 

1 

HQ  AFSC/DLCA 

2 

AMRL/THE 

3 

SAF/MXQ 

1 

SAMSO/DEV 

1 

AFIT/DE 

1 

HQ  AFSC/DEV 

1 

USN  Chief,  RSD/EQ 

1 

HQ  AFSC/SQPE 

1 

1 MSBH 

1 

HQ  ATC/SQPAP 

1 

6595  STESTG/TS 

1 

HQ  MAC/SOPE 

1 

USAF  Hoapital/SGPB 

1 

HQ  PACAF/SQPB 

1 

DOD  Exploaivaa  Safety  Board 

1 

HQ  EAC/SFQA 

1 

8AM80/WE 

1 

Cmdr,  UB  Army  Mad  Bioengrg  Lab 

1 

SAMTEC/WE 

1 

OEHL/CC 

3 

OEHL 

1 

AFWL/SUL 

1 

AFWL/WE 

1 

AFQL/XOP 

1 

HQ  AF8C/SD 

1 

L'SAFSAM/EDE 

2 

AFIT/Library 

1 

AFRPL/Library 

1 

Fadaral  Laboratory  Program 

1 

SAMTEC/SEH 

1 

USA  Chief,  RSD/EQ 

1 

OEHL 

1 

HQ  AU/AUL 

1 

SAMSO/SQ 

1 

HQ  USAFA/Library 

1 

AMD/RDU 

1 

HQ  AFTEC/SQ 

1 

ADTC/CSV 

1 

HQ  AFSC/SOSM 

■..I  . 

AFETR/DER 

1 

EPA/ORD 

1 

OASD/ (MRAGL) ES 

1 

HQ  NASA  (MHS-7) 

1 

Det  1 ADTC/CC 

1 

HQ  AFISC/SBS 

2 

AFCBC/SU 

1 

USAF  Hoapital/SGP , Eglin  AFB 

1 

Dot  1 ADTC/WE 

2 

1 STRAD/SEM 

1 

AFCEC/D1V 

1 

USAF  Hoepital/SGP , Edwards  AFB 

1 

Dat  1 ADTC/ECA 

5 

6395  STESTG/SZ 

1 

Dat  1 ADTC/ECW 

3 

SAX 

1 

AFTAC/TRA 

1 

Det  1 ADTC/EC 

1 

Dir/cc  Cham  Sya  Lab 

3 

OG-ALC/8GP 

1 

UBAF  Hoepital/SGPM 

1 

DDC/TCA 

2 

ARPA 

1 

Dafanaa  Raoh  a Bngrg/AD 

1 

USA  Bnvmtl  Hygn  Aganoy 

1 
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